Thermal entanglement of qubit pairs on the Shastry- Sutherland lattice 
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We show that temperature and magnetic field properties of the entanglement between spins on the two- 
dimensional Shastry-Sutherland lattice can be qualitatively described by analytical results for a qubit tetramer. 
Exact diagonalization of clusters with up to 20 sites reveals that the regime of fully entangled neighboring pairs 
coincides with the regime of finite spin gap in the spectrum. Additionally, the results for the regime of vanishing 
spin gap are discussed and related to the Heisenberg limit of the model. 
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I. INTRODUCTION 

In any physical system with subsystems in interaction, indi- 
vidual parts of the system are to some extent entangled, even 
if they are far apart, as realized already at the beginning of 
modern quantum mechanics sixty years ago. Today it has be- 
come appreciated that the ability to establish entanglement be- 
tween quantum particles in a controlled manner is a crucial in- 
gredient of any quantum information processing system 1 . On 
the other hand, it turned out that the analysis of appropriately 
quantified entanglement between parts of the system can also 
be a very useful tool in the study of many body phenomena, 
as is, e.g., the behavior of correlated systems in the vicinity of 
crossovers between various regimes or even points of quantum 
phase transition 2 . 

Quantum entanglement of two distinguishable particles in a 
pure state can be quantified through von Neuman entropy 3 ' 4 ' 5 . 
Entanglement between two spin-i particles - qubit pair - can 
be considered a physical resource, an essential ingredient of 
algorithms suitable for quantum computation. For a pair of 
subsystems A and B, each occupied by a single electron, an 
appropriate entanglement measure is the entanglement of for- 
mation, which can be quantified from the Wootters formula 6 . 
In general, electron-qubits have the potential for even richer 
variety of entanglement measure choices due to both their 
charge and spin degrees of freedom. When entanglement is 
quantified in systems of indistinguishable particles, the mea- 
sure must account for the effect of exchange and it must ade- 
quately deal with multiple occupancy states 7 ' 8 ' 9 ' 10 ' 11 ! 12 . A typ- 
ical example is the analysis of entanglement in lattice fermion 
models (the Hubbard model, e.g.) where double occupancy 
plays an essential roleii. 

In realistic hardware designed for quantum information 
processing, several criteria for qubits must be fulfilled^: the 
existence of multiple identifiable qubits, the ability to initial- 
ize and manipulate qubits, small decoherence, and the ability 
to measure qubits, i.e., to determine the outcome of compu- 
tation. It seems that among several proposals for experimen- 
tal realizations of such quantum information processing sys- 
tems the criteria for scalable qubits can be met in solid state 
structures consisting of coupled quantum dots 14 ' 15 . Due to the 



ability to precisely control the number of electrons in such 
structures 16 , the entanglement has become experimentally ac- 
cessible quantity. In particular, recent experiments on semi- 
conductor double quantum dot devices have shown the evi- 
dence of spin entangled states in GaAs based heterostuctures 17 
and it was shown that vertical-lateral double quantum dots 
may be useful for achieving two-electron spin entanglement 18 . 
It was also demonstrated recently that in double quantum dot 
systems coherent qubit manipulation and projective readout is 
possible^. 

Qubit pairs to be used for quantum information processing 
must be to a high degree isolated from their environment, oth- 
erwise small decoherence requirement from the DiVincenzo's 
checklist can not be fulfilled. The entanglement, e.g., between 
two antiferromagnetically coupled spins in contact with ther- 
mal bath, is decreased at elevated temperatures and external 
magnetic fiel d 20,21 ' 22 , and will inevitably vanish at some fi- 
nite temperature 23 . Entanglement of a pair of electrons that 
are confined in a double quantum dot is collapsed due to the 
Kondo effect at low temperatures and for a very weak tunnel- 
ing to the leads. At temperatures below the Kondo temper- 
ature a spin-singlet state is formed between a confined elec- 
tron and conduction electrons in the leads— For other open 
systems there are many possible sources of decoherence or 
phase-breaking, for example coupling to phonon degrees of 
freedom^ 5 .. 

The main purpose of the present paper is to analyze the ro- 
bustness of the entanglement of spin qubit pairs in a planar 
lattice of spins (qubits) with respect to frustration in magnetic 
couplings, elevated temperatures as well as due to increasing 
external magnetic field. The paper is organized as follows. 
Sec. |n] introduces the model for two coupled qubit pairs - 
qubit tetramer - and presents exact results for temperature and 
magnetic field dependence of the entanglement between near- 
est and next-nearest-neighboring spins in a tetrahedron topol- 
ogy. In Sec. [Ill] the model is extended to infinite lattice of 
qubit pairs described by the Shastry-Sutherland model 26 . This 
model is convenient firstly, because of the existence of sta- 
ble spin-singlet pairs in the ground state in the limit of weak 
coupling between the qubit pairs, and secondly, due to a rel- 
atively good understanding of the physics of the model in the 



thermodynamic limit. Entanglement properties of the Shastry- 
Sutherland model were so far not considered quantitatively. 
Neverteless, several results concerning the role of entangle- 
ment at a phase transition in other low-dimensional spin lattice 
system s 2 ' 27 ' 28 ' 29 ' 30,31 ' 32 , as well as in fermionic system s 33 ' 34,35 
have been reported recently. Near a quantum phase transi- 
tion in some cases entanglement even proves to be more ef- 
ficient precursor of the transition compared to standard spin- 
spin correlations 35 ' 36 . In Sec. |lV]we discuss entanglement be- 
tween nearest neighbors in the Heisenberg model, represent- 
ing a limiting case of the Shastry-Sutherland model. Results 
are summarized in Sec. [V]and some technical details are given 
in Appendix A. 



II. THERMAL ENTANGLEMENT OF A QUBIT 
TETRAMER IN MAGNETIC FIELD 

Consider first a double quantum dot composed of two adja- 
cent quantum dots weakly coupled via a controllable electron- 
hopping integral. By adjusting a global back-gate voltage, 
precisely two electrons can be confined to the dots. The inter- 
dot tunneling matrix element t determines the effective anti- 
ferromagnetic (AFM) superexchange interaction J ~ At 2 /U, 
where U is the scale of Coulomb interaction between two 
electrons confined on the same dot. There are several possi- 
ble configurations of coupling between such double quantum 
dots. One of the simplest specific designs is shown schemat- 
ically in Fig. Q2a): four qubits at vertices of a tetrahedron. In 
addition to the coupling A-B, by appropriate arrangements of 
gate electrodes the tunneling between A-C and A-D can as 
well be switched on. 

We consider here the case where J/U <C 1, thus double oc- 
cupancy of individual dot is negligible and appropriate Hilbert 
space is spanned by two dimers (qubit pairs): spins at sites 
A-B and C-D are coupled by effective AFM Heisenberg mag- 
netic exchange J and at sites A-C, B-C, A-D, B-D by J'. The 
corresponding hamiltonian of such a pair of dimers is given as 

Hi = J(Sa-S b +Sc-S d )+ (1) 
+ 2 J'(Sa • Sc + Sb • Sc + Sa • Sd + Sb • So) — 

- B(S Z A + S| + Sc + Sd), 

where Si = \gi is spin operator corresponding to the site i 
and B is external homogeneous magnetic field in the direction 
of the z-axis. Factor 2 in Eq. (HJ is introduced for convenience 
- such a parameterization represents the simplest case of finite 
Shastry-Sutherland lattice with periodic boundary conditions 
studied in Sec. [HI] 

A. Concurrence 

We focus here on the entanglement properties of two cou- 
pled qubit dimers. The entanglement of a pair of spin 
qubits A and B may be defined through concurrence 3 , C — 

2|cK-|-pQ!j_|_ — a||a^|, if the system is in a pure state |^ab) = 
Tl S s' a ss' \s)a\s')b, where \s)i corresponds to the basis | f )i, 
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Figure 1 : (Color online) (a) Two coupled qubit pairs (dimers) in tetra- 
hedral topology, (b) Shastry-Sutherland lattice as realized, e.g., in the 
SrCu2(B03)2 compound. 

| Concurrence varies from C = for an unentan- 

gled state (for example | f )aJ T)b) to C = 1 for com- 
pletely entangled Bell states 3 - ^(| T >a| T >b ± I 4 >a| 4 >b) 

or^(|T> A |4>B±|4>A|T>B). 

For finite inter-pair coupling J' ^ or at elevated tempera- 
tures the A-B pair can not be described by a pure state. In the 
case of mixed states describing the subsystem A-B the concur- 
rence may be calculated from the reduced density matrix /?ab 
given in the standard basis \ s)i\s')j 6 . Concurrence can be fur- 
ther expressed in terms of spin-spin correlation functions 2 ^, 
where for systems that are axially symmetric in the spin space 
the concurrence may conveniently be given in a simple closed 
form 37 , which for the thermal equilibrium case simplifies fur- 
ther, 

Cab = 2max(0, |(S+S B )| - yfipjpj) (p[p[)). (2) 

Here Sf' = (S^Y = Sf + iSf is the spin raising operator 
for dot i and P} = ±(1 + 25f), P/ = |(1 - 2Sf ) are the 
projection operators onto the state 1 1 )« or I 4 )i> respectively. 
We consider the concurrence at fixed temperature, therefore 
the expectation values in the concurrence formula Eq. (f2]) are 
evaluated as 

(0) = i^(n|0|n)e-^, (3) 

n 

where Z — J2 n e~ /3En is the partition function, (3 = 1/T, and 
{|n)} is a complete set of states of the system. Note that due 
to the equilibrium and symmetries of the system, several spin- 
spin correlation functions vanish, (S^S^) — 0, for example. 
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Figure 2: (Color online) (a) Zero-temperature concurrence Cab as a function of J'/ J and B/J. Different regimes are characterized by 
particular ground state functions \(f> n ) defined in Appendix A. (b) T/J — 0.1 results for Cab- (c) Next nearest concurrence Cac for T = 0, 
and (d) for T / J — 0.1. Dashed lines separate Cab(C) > from Cab(C) — 0. 



In vanishing magnetic field, where the SU(2) symmetry is 
restored, the concurrence formula Eq. (0 simplifies further 
and is completely determined by only one 38 spin invariant 



Cab = max(0, -2(S A ■ S B ) 



(4) 



The concurrence may be expected to be significant whenever 
enhanced spin-spin correlations indicate A-B singlet forma- 
tion. 



B. Analytical results 

There are several known results related to the model Eq. (Q}. 
In the special case of J' = 0, for example, the tetramer con- 
sists of two decoupled spin dimers with concurrence Cab 
(or the corresponding thermal entanglement) as derived in 
Refs. l2Cj||2ll Entanglement of a qubit pair described by the 
related XXZ Heisenberg model with Dzyaloshinskii-Moriya 
anisotropic interaction can be also obtained analytically 22 . 
Hamiltonian H4 with additional four-spin exchange interac- 
tion but in the absence of magnetic field was considered re- 
cently in the various limiting cases 39 . 

Tetramer model Eq. (Q~|i considered here is exactly solvable 
and in Appendix A we present the corresponding eigenvec- 
tors and eigenenergies. The concurrence C A b is for this case 



determined from Eq. (O with 

<S+S B ) = I [- eW/72 - CV 2 (e b + e- b ) /2 

+ e-'/W/e + e -''/2+tf (e b + e" 6 ) /4 
+ e--?/ 2 - 2 .?" (e 6 /4 + 1/3 + e~ b /4)] , (5) 

where j = j3J, f = j3J', b = /3B, and with 



/ pTl pTi\ _ 

\ A ^B / — ^ 



e-''/ 2 (e ±b ) + e- j / 2+Aj '/3 



/2 



Here 



+ e^' 2 -^' (1/6 + e ±b /2 + e ±2b ) 



Z = e 3j / 2 +2e j / 2 (e b +l + e- b ) 

+ e-j/2+4/ + e -j/2+2f (J + 1 + e -6J 

+ e -J/^' ( e 2b + e b +1 + e -b + e -2fe-) 



(6) 



(7) 



is the partition function. 

Alternatively, one can define and analyze also the entangle- 
ment between spins at sites A and C and the corresponding 
concurrence Cac can be expressed from Eq. © by applying 
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Figure 3: (Color online) (a) Temperature and magnetic field depen- 
dence of Cab for J'/ J = 0.4 and (b) J' — J. Dashed lines separate 
Cab > from Cab = 0. 



additional correlators with replaced B^C, 



(^A^C/ 



-j/2+4j' 



/3 



+ e - j ^- 2j ' (e fc /4 + 1/3 + e- b /4) 



(8) 



and 



C. Examples 

In the low temperature limit the concurrence is determined 
by the ground state properties while transitions between var- 
ious regimes are determined solely by crossings of eigenen- 
ergies, which depend on two parameters (J'/ J, B/J). There 
are 5 distinct regimes for Cab shown in Fig. Ufa): (i) com- 
pletely entangled dimers (singlets A-B and C-D, state \(p±) 
from Appendix A), Cab = 1; (ii) for B > J and smaller 
J'/ J the concurrence is zero because the energy of the state 
consisted of a product of fully polarized A-B and C-D triplets, 
|0i2), is the lowest energy in this regime; (iii) concurrence 
is zero also for J' > J/2 and low B/J, with the ground 
state |02). There are two regimes corresponding to | step in 
Cab where the ground state is either (iv) any linear combina- 
tion of degenerate states |<£e,7)> simultaneous A-B singlet 
(triplet) and C-D triplet (singlet) for J' < J/2, or (v) state 
\4>5) at J' > J/2 and larger B. Qubits A-C are due to special 
topology never fully entangled, and the corresponding Cac is 
presented in Fig. EJc). In the limit of J' 3> J the tetramer 
corresponds to a Heisenberg model ring consisted of 4 spins 
and in this case qubit A is due to tetramer symmetry equally 
entangled to both neighbors (C and D), thus Cac = \- 

At elevated temperatures the concurrence is smeared out as 
shown in Figs. |2fb,d). Note the dip separating the two dif- 
ferent regimes with Cab = 5, seen also in the Cac = 5 
case. This dip clearly separates different regimes discussed in 
the previous T = limit and signals a proximity of a disen- 
tangled excited state. For sufficiently high temperatures van- 
ishing concurrence is expected 23 . The critical temperature T c 
denoted by a dashed line is set by the magnetic exchange scale 
J, since at higher temperatures local singlets are broken irre- 
spectively of the magnetic field. 

A rather unexpected result is shown in Fig. (3ja) where at 
B > 2 J and low temperatures the concurrence slightly in- 
creases with increasing temperature due to the contribution of 
excited A-B singlet components that are absent in the ground 
state. Similar behavior is found for J' ~ around B ~ J, 
which is equivalent to the case of a single qubit dirtier^! 2 ! 
(not shown here). There is no distinctive feature in tempera- 
ture and magnetic field dependence of Cab when J' > J/2 
and a typical results is shown in Fig.[3jb) for J' > J. 



(P^PTI) 



e 3j / 2 /4 + e j/2 (l/2 + e ±b ) 
+ e -i/2+4i' /12 + e -i/ 2+2 i'e ±b /2 



-i/2-2/ 



(1/6 



72- 



(9) 



III. PLANAR ARRAY OF QUBIT PAIRS: THE 
SHASTRY-SUTHERLAND LATTICE 



Preliminaries 



The line 2 J' = J represents a particularly interesting spe- 
cial case where two dimers are coupled symmetrically form- 
ing a regular tetrahedron. An important property of this sys- 
tem is the (geometrical) frustration of, e.g., qubits C-A-B. 
Such a frustration is the driving force of the quantum phase 
transition found in the Shastry-Sutherland model and is the 
reason for similarity of the results for two coupled dimers and 
a large planar lattice studied in the next Section. 



The central point of this paper is the analysis of pair en- 
tanglement for the case of a larger number of coupled qubit 
pairs. In the following it will be shown that the results corre- 
sponding to tetramers considered in the previous Section can 
be very helpful for better understanding pair-entanglement of 
N > 4 qubits. There are several possible generalizations of 
coupled dimers and one of the simplest in two dimensions is 
the Shastry-Sutherland lattice shown in Fig. |TJb). Neighbor- 
ing sites A-B are connected with exchange interaction J and 
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next-neighbors with J'. The corresponding hamiltonian for 
N/2 dimers (N sites) is given with 

N 

H N = j^2Si- s, + j' s > ■ s j - B Yl s >- (10) 

{AB} {AC} i=l 

Periodic boundary conditions are used. For the special case 
N = 4 the model reduces to Eq. (Q~|i where due to periodic 
boundary conditions sites A-C (and other equivalent pairs) are 
doubly connected, therefore a factor of 2 in Eq. ([TJ, as men- 
tioned in Sec. HI1 

The Shastry-Sutherland model (SSM) was initially pro- 
posed as a toy model possessing an exact dimerized eigenstate 
known as a valence bond crystal 26 . Recently, the model has 
experienced a sudden revival of interest by the discovery of 
the two-dimensional spin-liquid compound SrCu2(B03) 2 4Q ' 41 
since it is believed that magnetic properties of this compound 
are reasonably well described by the SSM^ 2 -. In fact, several 
generalizations of the SSM have been introduced to account 
better for recent high-resolution measurements revealing the 
magnetic fine structure of SrCu2 (BOs) ^ 42 ' 43 ' 44 ' 45 . Soon af- 
ter the discovery of the SrCu2(BC»3)2 system, the SSM thus 
became a focal point of theoretical investigations in the field 
of frustrated AFM spin systems, particularly low-dimensional 
quantum spin systems where quantum fluctuations lead to 
magnetically disordered ground states (spin liquids) with a 
spin gap in the excitation spectrum. 

The SSM is a two-dimensional frustrated antiferromagnet 
with a unique spin-rotation invariant exchange topology that 
leads in the limit J ^> J' to an exact gapped dimerized ground 
state with localized spin singlets on the dimer bonds (dimer 
phase). In the opposite limit, J -C J', the model becomes 
ordinary AFM Heisenberg model with a long-range Neel or- 
der and a gapless spectrum (Neel phase). While two of the 
phases are known, there are still open questions regarding the 
existence and the nature of the intermediate phases. Several 
possible scenarios have been proposed, e.g.: either a direct 
transition between the two states occurs at the quantum criti- 
cal point near J'/ J ~ O.T 4 ^ 4 ^, or a transition via an interme- 
diate phase that exists somewhere in the range of J'/ J > 0.6 
and J'/ J < 0.9^. Although different theoretical approaches 
have been applied, a true nature of the intermediate phase (if 
any) has still not been settled. As will be evident later on, our 
exact-diagonalization results support the first scenario. 

The SSM phase diagram reveals interesting behavior also 
for varying external magnetic field. In particular, experiments 
on SrCu2(B03)2 in strong magnetic fields show formation of 
magnetization plateau s 41,49 , which are believed to be a con- 
sequence of repulsive interaction between almost localized 
spin triplets. Several theoretical approaches support the idea 
that most of these plateaus are readily explained within the 
(bare) SS M 46,50 ' 51 . Recent variational treatment based on en- 
tangled spin pairs revealed new insight into various phases of 
the SSM 48 . 

Although extensively studied, the zero-temperature phase 
diagram of the SSM remains elusive. This lack of reliable 
solutions is even more pronounced when considering thermal 
fluctuations in SSM as only few methods allow for the inclu- 
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Figure 4: (Color online) Results for the Shastry-Sutherland lattice 
with N = 20 sites and periodic boundary conditions. Presented are 
renormalized spin-spin correlation functions — 2(Sa • Sb,c) — § as 
a function of ,/'/ J and for various temperatures. Asterisk indicates 
critical J' c which roughly separates the dimer and Neel phase. 

sion of finite temperatures in frustrated spin systems. In this 
respect, the calculation of thermal entanglement between the 
spin pairs would also provide a new insight into the complex- 
ity of the SSM. 

B. Numerical method 

We use the low-temperature Lanczos method 52 (LTLM), an 
extension of the finite-temperature Lanczos method 53 (FTLM) 
for the calculation of static correlation functions at low tem- 
peratures. Both methods are nonperturbative, based on the 
Lanczos procedure of exact diagonalization and random sam- 
pling over different initial wave functions. A main advan- 
tage of LTLM is that it accurately connects zero- and finite- 
temperature regimes with rather small numerical effort in 
comparison to FTLM. On the other hand, while FTLM is lim- 
ited in reaching arbitrary low temperatures on finite systems, 
it proves to be computationally more efficient at higher tem- 
peratures. A combination of both methods therefore provides 
reliable results in a wide temperature regime with moderate 
computational effort. We note that FTLM was in the past suc- 
cessfully used in obtaining thermodynamic as well as dynamic 
properties of different models with correlated electrons as are: 
the t-J modelr 3 - the Hubbard modelr 4 . as well as the SSM 
models 

In comparison with the conventional Quantum Monte Carlo 
(QMC) methods LTLM possesses the following advantages: 
(i) it does not suffer from the minus-sign problem that usually 
hampers QMC calculations of many-electron as well as frus- 
trated spin systems, (ii) the method continuously connects the 
zero- and finite-temperature regimes, (iii) it incorporates as 
well as takes the advantage of the symmetries of the prob- 
lem, and (iv) it yields results of dynamic properties in the 
real time in contrast to QMC calculations where imaginary- 
time Green's function is obtained. The LTLM (FTLM) is on 
the other hand limited to small lattices which usually leads 
to sizable finite-size effects. To account for these, we ap- 
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Figure 5: (Color online) (a) Zero-temperature concurrence Cab for a 20-site cluster for various J'/ J and B/J. Shaded area represents the 
regime of fully entangled dimers, Cab = 1. (b) The corresponding results for T / J — 0.1. (c) Next nearest concurrence Cac for T = 0, and 
(d) for T I J = 0.1. Note qualitative and even quantitative similarity with the tetramer results, Fig.[2] Dashed lines separate Cab(C) > from 
Cab(c) = 0. 



plied LTLM to different square lattices with N = 8, 16 and 
20 sites using periodic boundary conditions (we note that 
next-larger system, N = 32, was too large to be handled 
numerically). Another drawback of the LTLM (FTLM) is 
the difficulty of the Lanczos procedure to resolve degener- 
ate eigenstates that emerge also in the SSM. In practice, this 
manifests itself in severe statistical fluctuations of the cal- 
culated amplitude for T — ► since in this regime only a 
few (degenerate) eigenstates contribute to thermal average. 
The simplest way to overcome this is to take a larger num- 
ber of random samples R 1, which, however, requires 
a longer CPU time. We have, in this regard, also included 
a small portion of anisotropy in the SSM (in the form of 
the anisotropic interdimer Dzyaloshinskii-Moriya interaction 
DZ J2 { ac}( s aS V c ~ S y A S£), D*/J ~ 0.01), which slightly 
splits the doubly degenerate single-triplet levels. In this way, 
R ~ 30 per S z sector was enough for all calculated curves 
to converge within ~ 1% for T/J < 1. Here, the number 
of Lanczos iterations M = 100 was used along with the full 
reorthogonalization of Lanczos vectors at each step. 



C. Entanglement 

Entanglement in the absence of magnetic field is most 
prominently reflected in spin-spin correlation functions, e.g., 



(Sa -Sb) and (Sa -Sc)- In zero temperature limit due to quan- 
tum phase transition at J' c these correlations change sign. In 
Fig. |4] are presented renormalized spin-spin correlation func- 
tions (for positive values identical to concurrence) as a func- 
tion of J' j J: (i) Cab > in dimer phase and (ii) Cac > 
in the Neel phase. Critical J' c is indicated by asterisk. The re- 
sults for N = 16 are qualitatively and quantitatively similar to 
the N — 20 case presented here. At finite temperatures spin 
correlations are smeared out as shown in Fig.|4]for various T. 
Limiting Heisenberg case, J' — > oo, is discussed in more de- 
tail in the next Section. J' = case corresponds to the single 
dimer limit 21 and Sec. [TTJ 

Complete phase diagram of the SSM at T = but with fi- 
nite magnetic field can be classified in terms of concurrence 
instead of spin correlations. In Fig. |3a)CAB is presented as 
a function of (J' / J, B / J) as in the case of a single tetramer, 
Fig. 12 a). Presented results correspond to the = 20 case, 
while N — 16 system exhibits very similar structure (not 
shown here). N = 8 and N — 4 cases are qualitatively sim- 
ilar, the main difference being the value of critical J' c which 
increases with N. Remarkable similarity between all these 
cases can be interpreted by local physics in the regime of fi- 
nite spin gap, J' < J' c . Qubit pairs are there completely entan- 
gled, Cab = 1) an d Cab ~ h f° r magnetic field larger than 
the spin gap, but B < J + 2J'. For even larger B concurrence 
approaches zero, similar to the N = 4 case. Concurrence is 
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Figure 6: (Color online) Temperature and magnetic field dependence 
of Cab for J'/J = 0.4 and N = 20. Note the similarity with 
the corresponding tetramer results, Fig. [3ja). Dashed lines separate 
Cab > from Cab = 0. 



zero also for J' > J' c , except along the B ~ 4J' line where 
weak finite concurrence could be the finite size effect. Similar 
results are found also for N — 16,8 cases, and are most pro- 
nounced in the N = 4 case. At finite temperature the structure 
of concurrence is smeared out [Fig. [2b)] similar to Fig.|2b). 

Concurrence Cac corresponding to next-nearest neighbors 
is, complementary to Cab. increased in the Neel phase of the 
diagram, Fig. [5jc). The similarity with N = 4, Fig. |2jc) is 
somewhat surprising because in this regime long-range cor- 
relations corresponding to the gapless spectrum of AFM-like 
physics are expected to change also short range correlations. 
The only quantitative difference compared to N = 4 is the 
maximum value of Cac ~ 0.3 instead of 0.5 (beside the crit- 
ical value J' c discussed in the previous paragraph). Concur- 
rence is very small for B > J + 2 J'. At finite temperatures 
fine fluctuations in the concurrence structure are smeared out, 
Fig. Ed). 

Temperature and magnetic field dependence of Cab in the 
dimer phase is presented in Fig.[6]for fixed J'/J = 0.4. Sim- 
ilarity with the corresponding N = 4 tetramer case, Fig. [3ja), 
is astonishing and is again the consequence of local physics in 
the presence of a finite spin gap. Finite size effects (in com- 
parison with N — 16 and N — 8 cases) are very small (not 
shown). Dashed line represents the borderline of the Cab = 
region: critical T c as 0.75 J valid for Bj J < 3, that is in this 
regime nearly independent of B, is slightly larger than in the 
single tetramer case where its insensitivity to B is even more 
pronounced. 



IV. HEISENBERG LIMIT 

The concurrence corresponding to next-nearest neighbors 
in SSM, Cac. is non zero in the Neel phase for J' > J' c . 
Typical result for concurrence in this regime (for fixed J'/J = 
1) in terms of temperature and magnetic field is presented in 
Fig. ITJa). At zero temperature the concurrence is zero for B > 
4 J' [compare with Fig.[2|c) and Fig. [2c)]. 




Figure 7: (Color online) (a) Next nearest neighbor concurrence Cac 
for J' = J. (b) Heisenberg lattice result as a special case of the 
SSM, J = 0. Shaded region represents Cac = 0. In the line shaded 
region (low finite temperature and large magnetic field) our numer- 
ical results set only the upper limit Cac < 5 • 10~ 4 . Dashed lines 
separate Cac > from Cac = 0. 

In the limit J = the model simplifies to the AFM Heisen- 
berg model on a square lattice of N sites, 

N 

{AC} i=l 

Several results for this model have already been presented for 
very small cluster s 55 ' 56,57 , however the temperature and mag- 
netic field dependence of the concurrence for systems with 
sufficiently large number of states and approaching thermo- 
dynamic limit has not been presented so far. 

In Fig.[7lb) we further presented temperature and magnetic 
field dependence of concurrence for the Heisenberg model for 
N = 20 (results for N = 16 are quantitatively similar, but not 
shown here). Temperature and magnetic field dependence of 
Cac exhibits peculiar semi-island shape where at fixed value 
of B the concurrence increases with increasing temperature. 
This effect is to some extent seen in all cases and is the con- 
sequence of exciting local singlet states, which do not appear 
in the ground state. At T — ► finite steps with increasing B 
correspond to gradual transition from the singlet ground state 
to totally polarized state with total spin 5=10 and vanish- 
ing concurrence. This is in more detail presented in Fig.^a) 
for various N = 4, 8, 16, 20. At B = and for N = 20 
we get Cac = 0.19. It is interesting to compare this results 
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Figure 8: (Color online) (a) Zero-temperature concurrence Cac in 
the Heisenberg limit as a function of B/J' and for various N — 4, 
8, 16, 20. Sections with different total spin values are addition- 
ally labeled, (b) Finite-size scaling of concurrence in the absence of 
magnetic field. Full line represents the fit corresponding to Ref . I58L 
C AC ~| + 2iV- 3 / 2 . 



with the known finite-size analysis scaling for the ground state 
energy of the Heisenberg model^. The same scaling gives 
Cac ~ | + 27V -3 / 2 . Our finite-size scaling, Fig. Ob), is in 
perfect agreement with this result for N — ► oo at T — and 
B = 0. 

In the opposite limit of high magnetic fields, the vanishing 
concurrence Cac = 0, is observed for B above the critical 
value B c = 4 J' for all system sizes shown in Fig.[8|a). This 
result can be deduced also analytically. Since in a fully po- 
larized state Cac = 0, this B c actually denotes a transition 
from Si = N/2 — 1 to Sq = N/2 ferromagnetic ground 
state with energy Eq = N(J' — B)/2. The energy of the 
one-magnon excitation above the ferromagnetic ground state 
is given by the spin wave theory, which is in this case exact, as 
Ei = Eq — J'(2 — cos k x a — cos k y a) + B, where (k x , k y ) is 
the magnon wave vector and a denotes the lattice spacing. Ev- 
idently, a transition to a fully polarized state occurs precisely 
at B c = A J' at (it /a, ir/a) point in the one-magnon Brillouin 
zone. 



V. SUMMARY 

The aim of this paper was to analyze and understand 
how concurrence (and related entanglement) of qubit pairs 
(dimers) is affected by their mutual magnetic interactions. In 
particular, we were interested in a planar array of qubit dimers 
described by the Shastry-Sutherland model. This model is 
suitable due to very robust ground state composed of entan- 
gled qubit pairs which breaks down by increasing the in- 
terdimer coupling. It is interesting to study both, the en- 
tanglement between nearest and between next-nearest spins 
(qubits) at finite temperature and magnetic field. The results 
are based on numerical calculations using low-temperature 



Lanczos methods on lattices of 4, 8, 16 and 20 sites with pe- 
riodic boundary conditions. 

A comprehensive analysis of concurrence for various pa- 
rameters revealed two general conclusions: 

(1) For a weak coupling between qubit dimers, J' < J' c , 
qubit pairs are locally entangled in accordance with the local 
nature of the dimer phase. This is due to a finite singlet-triplet 
gap (spin gap) in the excitation spectrum that is a consequence 
of strong geometrical frustration in magnetic couplings. The 
regime of fully entangled neighbors perfectly coincides with 
the regime of finite spin gap as presented in Fig. [9] Calcu- 
lated lines for various system sizes N in Fig. [9ja) denote re- 
gions (shaded for A" = 20) in the (J 7 /J, B/J) plane where 
Cab = 1 at T = 0. In the lower panel [Fig. fyb)] the lines 
represent the energy gap E\ ~ Eqs between the first excited 
state with energy E\ with total spin projection S z — 1 and 
the ground state with energy Eqs and total spin projection 
S z = 0, calculated for B = 0. For J' < J' c (full fines) 
Ei — Eqs corresponds to the value of the spin gap. With an 
increasing magnetic field the spin gap closes (shaded region 
for N — 20) and eventually vanishes at the Cab = 1 border 
line. Shaded regions in Figs. |9|a),(b) therefore coincide. Note 
also that the results for A" = 16 and 20 sites differ mainly in 

J'c- 

As a consequence of finite spin gap and local character of 
correlations it is an interesting observation that even N = 4 
results as a function of temperature and magnetic field quali- 
tatively correctly reproduce A^ = 20 results in the regime of 
J' < J' c . The main quantitative difference is in a renormal- 
ized value of J' c = J/2 for N = 4, as is evident from the 
comparison of Figs l2l3l and Figs l5l6l This similarity of the re- 
sults appears very useful due to the fact that concurrence for 
tetrahedron-like systems (N — 4) is given analytically (Sec. 

hd. 

(2) In the opposite, strong interdimer coupling regime, J' > 
J' c , the excitation spectrum is gapless and the concurrence be- 
tween next-nearest qubits, Cac> exhibits a similar behavior as 
in the antiferromagnetic Heisenberg model J/ J' — > 0. Our 
B = results coincide with the known result extrapolated to 
the thermodynamic limit Cac ~ In finite magnetic field 
and T = the concurrence vanishes at B c = 4 J' when the 
system becomes fully polarized (ground state with the total 
spin S = N/2). However, at elevated temperatures the con- 
currence increases due to excited singlet states and eventually 
drops to zero at temperatures above T c w J'. 

We can conclude with the observation that our analysis 
of concurrence and related entanglement between qubit pairs 
was also found to be a very useful measure for classify- 
ing various phases of the Shastry-Sutherland model. As our 
numerical method is based on relatively small clusters, we 
were unable to unambiguously determine possible interme- 
diate phases of the model in the regime J' ~ J' c , but we be- 
lieve that concurrence will prove to be a useful probe for the 
classification of various phases also in this regime using alter- 
native approaches. However, we were able to sweep through 
all other dominant regimes of the parameters including finite 
temperature and magnetic field. 



1 

0. 

> 0.6 
) 

0.4 

0.2 


1 

5 0.8 
IS 0.6 

""" 0.2 
0. 




S = 1, s z 



' 1 1 1 1 .• 4 1 


1 I 1 

16 


; \ \ 

singlet - triplet gap 

" (b) 

i.i.i 


• (V=20 

1 



0.2 



0.4 0.6 

ru 



0.8 



Figure 9: (Color online) (a) Zero-temperature Cab = 1 region in the 
plane (J' / J, B/J) for various N. (b) The corresponding spin gap 
at B — (the energy of the lowest total S z — 1 state relative to the 
ground state energy). 
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Appendix A: EIGENENERGIES AND EIGENVECTORS FOR 
PERIODICALLY COUPLED TWO QUBIT DIMERS 

Consider two qubit dimers coupled into a tetramer and de- 
scribed with the Hamiltonian Eq. (Q]i and Fig.QJa). The model 
is exactly solvable in the separate {S, S z } subspaces corre- 
sponding to different values of the total spin S and its z com- 
ponent S z . Following the abbreviations for singlet and triplet 
states on nearest-neighbor (dimer) sites i and j, 

1 



14) 



= I Vilj + litj)) 

= ITiTi), 

= | lilj), 



(Al) 



the resulting eigenstates |0fc) and eigenenergies Ek corre- 
sponding to the hamiltonian Eq. (HJ are: 



S = Q,S' = 0: 

\<t>i) = |sab)|scd), 
E x = -3J/2, 

|02> = ^(-|*AB>I*&D) 



rAB/ rCD/ 



Eo 



J/2 - 47'. 



(A2) 

l*AB>l*CD» 

(A3) 



103) = |sAb)|* cd ), 

104) = |£ab)I s cd), 
E 3A = -J/2 - B, 




(A4) 


i&> = ^(i*ab>i*ct>-i*Ib 


>I*&D»> 




E 5 = J/2 - 2 J' - B. 




(A5) 


S = 1,S Z = 0: 






|0 6 ) = |«ab)|*cd>» 

107) = |*AB>I*CD>, 

E 6 . 7 = -J/2, 




(A6) 


I08) = 4=(|4b)IW-|*Ib 

V2 


)I4d))> 




E$ = J/2 - 2 J'. 




(A7) 


S = 1, S z = 1 : 






109) = |sab)|*c d ), 






I0io) = I<ab)I s cd), 
-^9,10 = — J/% + B, 




(A8) 


|0h) = ^(-|4b)I4d) + I* 


ab)I*cd))> 




En = J/2 - 2 J' + B. 




(A9) 


S — 2, S z = -2 : 






|0i2) = |*ab)I*cd)« 
E12 = J/2 + 2 J' - 2B. 




(A10) 


S = 2,S Z = -1: 






I0i 3 ) = ^(|*ab)I*cd) + I*ab 


>I*CD», 




E 13 = J/2 + 2 J' - B. 




(All) 


S = 2,S Z =Q: 






1014) = |(2|^ B >I*CD> + I4b) 


l*CD) + l*AB> 




E U = J/2 + 2 J'. 




(A12) 


S = 2, S z = 1 : 






1015) = ^=(l*AB>fe> + l*AB 


)I*&d)), 




Ei 5 = J/2 + 2 J' + B. 




(A13) 



S = 2,S Z = 2 

^AB/l"( D 



1016) = |<AB>I*CD>> 
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Exe = J/2 + 2 J' + 2B. 



(A14) 
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